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AUTHOR’S PREFACE 


High school students, especially those interested in mathematics, 
physics and engineering, often ask, ‘What is “higher” mathematics? 
Sometimes they discuss this and similar questions at mathematics 
clubs at schools. 

In this book I have tried to explain, in a way a high school 
pupil would understand, certain concepts of higher mathematics *) 
such as the derivative, differential equation, the number e, and 
natural logarithm (pupils are more apt to be aware of and interested 
in the latter two concepts). Wherever possible, I have tried 
to illustrate the concepts with problems taken from _ physics. 
In addition, I have tried to show that the concepts of “higher 
mathematics” are mathematical reflections of actual processes, that 
mathematics and life are connected, not separated, and that 
mathematics is a growing, not an unchanging, completed science. 

Not all proofs and arguments are presented with complete mathe- 
matical rigour. Some arguments are presented for illustration. 
This method seems to me more appropriate for a general book. 

The book can be used by mathematics and physics clubs at 
school. Part of the material is taken from lectures the author 
gave at the request of the advisers of school mathematics clubs 
at the Moscow State University. 

I would like to take this opportunity to express my sincere 
gratitude to A. I. Markushevich and A. Z. Ryvkin for their valuable 
advice and comments on the manuscript. 


*) The reader may become acquainted with certain concepts of higher 
mathematics by reading other books of this series: A I. Markushevich 
Ploshchadi i logariphmy (Squares and Logarithms) (No. 9) and I. P. Natanson 
Summirovanie beskonechno malykh velichin (The Summation of Infinitesimals) 
(No 12) 
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The Problem of a Free Falling Body 


Statement of the Problem 


We shall first estimate the velocity of a body that falls to the 
ground from a certain height. 

Elementary physics tells us that a body falling freely in a vacuum 
obtains after ¢ seconds from the start of the process the velocity 


v=v9+gt (1) 


where vg is the initial velocity and g the gravitational acceleration. 

In certain cases formula (1) remains approximately valid for 
a body falling in the atmosphere (as distinct from a vacuum), 
in other cases it can lead to gross errors. For example, it is 
applicable to a body falling a short distance. When a body falls 
from a great height, the expression of velocity greatly differs from 
formula (1). In 1945 V.G. Romanyuk in a free-fall jump fell 
12000 m before he opened his parachute. In a vacuum a body 
which fell from such a height (without initial velocity) would at 
ground level have attained a speed of around 500 m/s. In fact, 


2 
it follows from the formula s = x that falling time equals (in 


/=~ 2x 12000 m 
= ~ 49, 
= xeon 9.8 m/s? mak 


and from (1) we obtain the speed 
v=gt ~9.8 m/s? x 49.5 s = 485 m/s 


(We could have used the formula v? = 2gs right away.) At the same 
time it has been established that the speed of the parachutist 
in the course of a free-fall jump reaches 50-60 m/s and does not 
exceed this value. Thus, in this case formula (1) leads us to an 
erroneous conclusion. 

Here is another example: the parachute is designed so that a man 
descending with it approaches the ground with a speed of 


a vacuum) 
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approximately 6.5 m/s no matter from what altitude the descent 
is being made. 

Clearly, formula (1) is in this case of no value. 

This points to the conclusion that the speed of a body falling 
in the atmosphere reaches a constant value in the course of time. 
In other words, some time after the start of the fall the motion 
of the body becomes uniform and its acceleration becomes zero. 
This means that the resultant force (the sum of all forces) acting 
on the body is zero. 

It is easily seen why formula (1) is of no use for calculating 
the velocity of a body falling in the atmosphere. The formula is 
deduced on the assumption that the body moves under the action 
of only one force, namely, the force of gravity 


P=mg (2) 

On the other hand, we have established that for a body falling 
in the atmosphere the resultant force (some time after the motion 
has begun) becomes zero, i. e. the force of gravity P is compensated 
by some other force which was not taken into account when formula 
(1) was deduced. This compensating force is the air resistance, 
or drag. It is the drag which prevents the parachutist from attaining 
excessive speed — it “supports” him, so to speak. 

How can drag be taken into account? Suppose there is no wind. 
For a stationary body the drag is zero. The faster the body moves 
the more difficult it becomes for it to penetrate the air, i.e. the 
drag increases. This fact may be easily demonstrated on a still day 
if you increase your speed: start walking, then running, cycling, etc. 
We shail suppose that the magnitude of this force is proportional 
to the velocity, i.e. that it equals bv, where b is the proportionality 
factor and v the speed of motion. This assumption is well supported 
by experiment for low speeds*) not exceeding 1-2 m/s. The factor b 
depends on the dimensions and on the shape of the body. 
For instance, the speeds being equal, the drag of a ball is 
approximately 20 times that of a “cigar-shaped” body of equal 
cross section (Fig. 1). 

Confining ourselves to these brief remarks, in future we shall 
assume that the drag (for which we shall adopt the designation S) is 

S = —bv (3) 
the minus sign shows that the drag has a direction opposite to that 
of velocity. 


*) Note that for speeds above 1-2 m/s the drag exceeds the value bv. 
It is sometimes assumed that the drag is proportional to the square 
of velocity 
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Fig 1 

All said, we shall assume that a body thrown downwards with 
an initial velocity is acted upon by only two forces: the force 
of gravity P and the drag S. Making use of Newton’s Second Law 
we may write: 


ma=P+S (4) 
where m is the mass of the body and a its acceleration. It is more 
convenient to point the positive direction vertically downwards than 
upwards since the velocity of a falling body is directed downwards, 
and, if this convention is adopted, it will be positive. The force 
of gravity which is also directed downwards will also be positive. 
On the other hand, the drag will have the direction opposite 
to the velocity, i.e. upwards, and will for this reason be negative. 
Hence, substituting the values (2) and (3) of P and S into formula 
(4) we obtain 


or 


a= -2(o- ™) (5) 


It is natural to assume the acceleration to be positive if it is 
directed downwards and negative if it is directed upwards. 

Equation (5) relates the acceleration to the speed of motion which 
we do not know yet. It is this equation which we shall have 
to use to determine how the velocity of the moving body will 
change in time. 


Qualitative Solution of the Problem 


As a result of the above discussion we have obtained equation 
(5) for the velocity of a falling body. Now it is up to us to solve 
the equation. For this reason the following arguments are of 
a purely mathematical nature, although to help the reader see 
the point we shall continue to treat the matter of the velocity 
of a falling body. 


Equation (5) connects two unknown quantities: velocity and 
acceleration. Assigning an arbitrary value to the acceleration we 
shall, using equation (5), obtain a definite value for the velocity. 
Therefore, our first impression is that one equation (5) is not enough 
to find two values v and a. 

However, this impression is erroneous. The acceleration of a body 
is completely determined by the manner in which its velocity 
changes with time. Therefore equation (5) contains not completely 
arbitrary values, but related values v and a. This makes it possible 
to solve equation (5). Studying the relation between the velocity 
and the acceleration, we shall later on arrive at the concept of 
the derivative. 

We shall prove two properties of velocity which follow from 
equation (5); these properties will make the process perfectly clear 
(in the assumptions made above). In due course we shall also 
obtain an accurate formula for the velocity. 

PROPERTY 1 If initially the velocity vo of a falling body is less than 


ae then in the course of motion the relation v< it always 


holds. If; on the other hand, vg > ae then v> re will always be 
true. 
mg 


Suppose this is not so. For instance, let v9 < ari and at a cer- 
tain instant ft, (i.e. t; seconds after the body started falling) let 
the velocity exceed o It follows that at some intermediate 
moment (perhaps, on numerous occasions) the velocity was equal to 
ae Let to be the last instant (in the first t; seconds) when the 
mg 

b 
for the interval tp —t,. It follows from formula (5) that during 
all this time the acceleration a was negative. However, this 
contradicts the statement that within the corresponding interval 
Ne. 

b 


proves that the velocity cannot exceed the value of 


velocity was equal to , so that the inequality v> ai held 


of time the velocity has exceeded . The resulting contradiction 


Bcd 

= 
The case vp > _ should be considered in the same way. 
a the velocity of a falling body increases 


with time, all the time approaching the value _ if, on the other 


PROPERTY 2 If vo < 
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hand, vo > aa the velocity decreases steadily also approaching the 


value “2 
Evidently, if, for instance, vo > a the relation v at will, 


as we know from Property 1. hold for the motion as a whole. 
It follows from formula (5) that the acceleration will be negative 
and therefore the velocity will steadily decrease. 


Let us prove that in the course of time the difference v — = 


will become less than any predetermined small quantity h (this 
may be, for instance, made equal to 0.001 m/s). Therefore, consider 


the instant 
(+ Lt), 
t*= b 


hb 
During the time from the start of the motion up to the instant 
t* the velocity of the falling body has decreased from the value vo 


to a value not less than ™, ie. it has decreased no less than by 


an amount ee 
Vo — mg = —t* 
b m 
It follows that at some intermediate instant the acceleration 
hb ; 
was not greater than —, for should the acceleration have been 
m 
greater during the whole interval of time the velocity would have 
hb 
decreased more than by the amount —t*. 
: ; hb 
So let us assume that at the instant t’ |a| << — 
m 


In accordance with (5) it follows that 
— | My gy gm hb 
aaa; | a ae 


i.e. at the instant t’ the difference between the velocity and the 
value a is less than h. This will hold for all future instants as 
well, for the velocity continues to decrease always remaining greater 


than oe: 


b 
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It should be noted that we have thus proved a somewhat more 
precise statement than Property 2, namely, that not later than 
after the time 


Fk pnt | 

t* =| b | hb S (6) 
after the fall has started the difference between the velocity and 
the value “Z will be less than h. 


Properties 1 and 2 provide us in a certain sense with a solution 
of the problem we have set ourselves. Although we have not yet 
obtained a precise formula for the velocity, we have, nevertheless, 
found out qualitative laws for the variation of the velocity, i.e. we 
know now how it is going to change with time. 

For instance, let us observe the motion of the parachutist. 
If he opens his parachute the moment he jumps, his velocity will 
increase from zero, but shall never exceed the value “A The 
value mg (the weight of the jumper with the parachute) is known, 
and b depends on the parachute’s diameter. This makes it possible 
to calculate the dimensions of the parachute so that the maximum 


possible speed of descent equal to a guarantees a safe landing 


for the parachutist. On the other hand, in case of a free-fall jump 
(with the parachute tucked away) the coefficient for the drag will 
have a new value, which we shall denote by b’ and which is 
smaller than in case of a descent with an open parachute For 
this reason the maximum speed a is greater in the latter case. 
It follows that before the parachute is opened in the course 


of a free-fall jump the parachutist’s speed will be greater than 
mg 


oe and, according to Property 1, after the parachute is opened 


it will decrease approaching _ at the same time remaining 


greater than a Thus in this case, too, some time after the 
parachute has been opened the jumper is able to make a safe 


landing. 

Let us take a numerical illustration. 

Example 1. Let a parachute be designed so that a jumper’s speed 
with the parachute open approaches the limiting value of 6 m/s, 
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Le. = 6 m/s. The jumper opens his parachute while falling 


at a speed of 50 m/s. How much time will lapse before his speed 


becomes less than 10 m/s, i.e. differs from the limiting value 
a =6 m/s by an amount less than h=4 m/s? 


Solution. From the equation — = 6 m/s we obtain: 


mim i Ooms _ 
b  ~=6b g~~—s«i10 m/s? Bie 
Next using formula (6) we conclude that the speed of descent 


will differ from the limiting value “a =6 m/s by h=4 m/s not 


; 1 : : 
later than after a time [ v9 - BES seconds, i.e. with the 


assumptions made, after 
1 


(50 m/s — 6 m/s) x 0.6 s x ami = 6.6 s 


Formula for the Velocity 
of a Falling Body. 
The Number e 


Properties 1 and 2 demonstrate how the velocity of a falling 
body changes with time. In this section we intend to obtain 
an exact formula for the velocity of a falling body. A certain 
number enters into the expression for the velocity whose value 
within the fifth decimal place is 2.71828... Frequently used in 
“higher” mathematics, this number is denoted by the letter e 
(in the same way as another frequently used number 3.14159..., 
which expresses the ratio of the circumference to the diameter, 
is denoted by the letter x). Later we shall see why this number 
e = 2.71828... enters into the formula for the velocity and how it 
should be accurately calculated. For the present, we shall write 
the formula for the velocity of a falling body without deduction 
and shall consider some examples illustrating its application. 

let v9 be the initial velocity of a falling body and v, the velocity 
of the body at the instant t (i.e. t seconds after the body has 
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started falling). Then we have 
4 
= (m-B)e m= 


This is the exact solution of equation (5). The proof of formula 
(7) will be given below. Let us consider some examples. 

Example 2. Let us demonstrate that formula (7) immediately 
leads to the qualitative laws for the variation of the velocity 
(Properties 1 and 2) obtained above. 


b 
The number e , which is obtained by raising the number e 
to a negative power, is, in effect, positive and less than unity, 


_» _4 _»y 
ie. 0<e ™<1. Ast increases, e ™ = (. 4) decreases (and 
may be made for sufficiently large t’s as small as desired). Therefore 


it follows from formula (7) that, for instance, for v9 > 2 the 


velocity v, always exceeds the value = (tr Vo — a > 0) and 


that it decreases with time approaching =. 


Example 3. Let us calculate, using formula (7), the velocity 
of the parachutist 6.6 seconds after he opens his parachute in 
a free-fall jump. We shall use the same numerical values as 
in Example 1 (p. 12), i. e. “A = 6 m/s, Vo = 50 m/s. (We have seen 
that this velocity must be less than 10 m/s.) 

Solution. We have 


Es Ie = ee 
Sa re = 10 m/s*:6 m/s = - 


Next, using the tables of logarithms (the common logarithm of 
the number e is approximately 0.4343), we easily obtain 


( -4,) b 5 
logio emjy=—- —, tlogioe y = oe x 6.6 x 0.4343 = 


= —4.7773 = —5 + 0.2227 = 5.2227 


whence 


b 
e ™' ~ 0,0000167 


Substituting this value into formula (7), we obtain 
66 ,~ 6 m/s + (50 m/s — 6 m/s) x 0.0000167 = 6.000735 m/s 


Similarly, with the aid of formula (7), it may be easily calculated 
that the parachutist speed will be 10 m/s (the conditions being 
3 logig 11 

5 logige 

Thus, when the parachute in a free-fall jump is opened, the speed 
of descent decreases in 1-2 s from 50-60 m/s almost to the normal 
speed of 6-7 m/s under an open parachute. Here the parachutist 
slows down considerably, i.e. undergoes the great force of the 
upward jerk of the parachute, caused chiefly by air resistance. One 
who has watched free-falls (at aviation holidays, for example) may 
have seen how a person rapidly plunging to the ground, suddenly, 
the moment the parachute opens, slows down sharply; it may even 
seem that for an instance he dangles in midair. 

Example 4. Let a parachutist’s rate of descent in a free-fall 


the same) ft = s = 1.44 s after the parachute opens. * 


jump approach the limiting value of a = 50 m/s. We shall set 


the speed at the start at zero. What will the error be if instead 
of formula (7) we use formula (1), applicable for a body falling 
in a vacuum? 

Solution We have 


=P = 8 10 m/s? 72 
mm mg 50 m/s s 
b 


Thus, in accordance with formula (7), the parachutist’s rate of 
descent will be 
v, = 50(1 — e ° 74) 


We obtain the following value of the speed of a body falling 
in a vacuum from formula (1). 


v, = gt = 10t 
Thus the ratio of the velocities becomes 
POUR Pipa 
gt t 


* The speed of descent will in fact approach the limiting value of 
6 m/s still more rapidly since expression (3) for the drag is really valid 
only for small velocities. For greater speeds the drag’s dependence 
on speed is greater than bv 
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Putting t= 1 s, we obtain (after some easy calculations with the 
aid of the logarithmic tables) for the above ratio the value ~ 0.91, 
and for t=2 s—the value ~ 0.82. We see that already the first 
few seconds after the body has started falling, its speed (because 
of the drag) differs noticeably from the value gt. 

Now let us prove formula (7). Therefore, let us first elucidate 
the relation between the velocity and the acceleration. If v, is the 
velocity of the body’s motion at the instant t, and v,,, its velocity 
h seconds after this instant (i.e. at the instant ++ h), the ratio 
rth — YY 

h 
interval h and denoted by a,,: 


is termed the average acceleration of the body in the 


Dr+h — Ur 

ay = h 
If h is very small (say, 0.01 s, or smaller, depending on the 
nature of the body’s motion) the acceleration changes little during 
such a small interval, so that a,, differs little from the value a, 
at the instant t. The difference between a, and a,, is the smaller 
the less is h. In other words, should we put h progressively 
smaller (say, 0.1 s, 0.01 s, 0.001 s, etc.), leaving t constant, a,, 
would change approaching a, closer with every step. In mathematics 
this fact is written down as follows: 
Vr+h — Ys 


a, = lima,, = lim 
h-0 h-O0 


The sign lim means the limit of the expression which follows 
it (i.e. the expression a,,); the designation h-—0 below indicates 
that we have in mind the limit a,, when h tends to zero. 

Thus we have obtained the relation expressing the dependence 
of the acceleration on the velocity. Let us now prove three 
additional properties of the velocity of motion we have been 
discussing. These properties will help us to prove formula (7). 

PROPERTY? 3 If the velocity and the acceleration of a moving body 
satisfy relation (5), this means that further variation of the velocity 
is uniquely determined by the initial velocity vo. 

Suppose the contrary is true. Let two bodies T and T* with 
equal values of m and b move in such a manner that their 
velocities and accelerations satisfy relation (5). Suppose at the 
instant t =0 these bodies had equal initial velocities vg, but tr, 
seconds later their velocities became different. say the velocity 
of the first body v, became greater than the velocity of the second 
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body vf. To be definite, suppose that vg > ~ (should the op- 


posite hold, the proof would be quite similar). Let to be the last 
instant (during the first t,; seconds) when the velocities of both 
bodies were equal. In that case in the period t, — to the velocity 
of the first body v was all the time greater than the velocity 
of the second body v*, i.e. v > v*. Hence it follows that 


8 Se 8 
v B27? ' 
both quantities v — = and v* — st being positive since 


v9 > 2 (see Property 1). In accordance with formula (5) we 
conclude from the inequalities 


8 ye 
v 5 >v b >0 


that the accelerations of both bodies a and a* are negative, the 
magnitude of a being greater than that of a*. This, however, means 
that the velocities of the bodies T and T* have decreased during 
the period t, — tf , the velocity of the body T having decreased 
more than that of the body T*, i.e. at the instant t, the velocity 
v should be less than v* (because at the instant ty these velocities 
coincided). We, however, made the opposite assumption. The result- 
ing contradiction proves that Property 3 is correct. 

PROPERTY 4. If two identical* bodies T and T* start falling 
simultaneously with initial velocities vg and v§, their velocities v, and 
v¥ at any instant t will satisfy the relation 


we — pg — 8. 


b beds ce SP (8) 
mg mg 
Be pe ee ae 


To prove this consider an imaginary body T which moves so 
that at the instant ¢ its velocity is equal to 


= 40, + (1 — 


In the sense that their respective m's and b’s are equal. 


.! 
~ 


where 


We shall demonstrate that the velocity and the acceleration of 
this imaginary body will satisfy relation (5). Find the average 
acceleration G,, of this imaginary body inside the interval t, t + h. 
We have 


~ gl |_ — gy 8. 
7 ee : ie + (1-4) | [asta q) “| 
h 


v — 
= qth = 4a,, 


where «i,, is the average acceleration of the body T inside the same 
interval. Should we take progressively smaller values of h in the 
relation d,, = qa,,, @,, would approach the acceleration 4, of the 
imaginary body at the instant t. and a,, the acceleration a, of the 
body T at the same instant. In this way we shall obtain 4, = qa, 
(for an arbitrary instant t) and with the aid of equation (5) 


i.e. relation (5) is satisfied for the imaginary motion we have been 
discussing above. 7 
Next the initial velocity of the imaginary body T is equal to 


Bo = 9 +1 — a) ME = (v9 ~ GE) + = 
iar 
75 (mt), Mo 
mg b b 
vg — — 


b 
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Thus both bodies T and T* have identical initial velocities and 
both move so that their velocities and accelerations satisfy equa- 
tion (5). In accordance with Property 3 it follows that their 
velocities 3, and v*¥ coincide at every instant ft, i.e. 


uf = qu, + (1 — ge 


Hence we obtain 


t-B [msa-o BR] oe F) 


= =q= 


a8 


thereby proving Property 4. 


PROPERTY 5 For any instants t and t the relation holds 


a (9) 


where U9, Uys Up Up4, are the velocities of the falling body T at the 
instants 0, t, t, t +1. 

Indeed, start observing how the body Tis falling from the instant 
t onwards. In t seconds (i.e. t+ seconds after the motion has 
started) the velocity will be v,.,. This means that if at the instant 
t =0 we in addition to the body T drop another body T* whose 
initial velocity v} equals v,, the velocity of the latter body v¥ at 
the instant ¢ will be equal to v,,,, i.e. v¥ = v,4,. Hence we obtain 
from (8) 


mg mg 
Pree Ee 
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or otherwise 


(ou B)(o-)> (6-8) 6-8) 


Dividing both parts of the equation thus obtained by 
2 
vo — = , we obtain the required relation (9). 
After deriving formula (9) let us start to calculate the exact 
value of the velocity v,. In order to avoid lengthy formulas let us, 


as a temporary measure, introduce the notation 


= 
u, = 
go eee 
ob 
thus simplifying formula (9): 
U4, = UU, (10) 


For t =¢t formula (10) yields 
uz, = ur 

Similarly, putting t = 2r, we obtain from (10) 

Us, = Ulla, = Uy = UP 
and for t = 3t the result will be 

Ug, = U3, = uu? = ue 
etc. Continuing in the same way, we shall find that for any 
positive integer n the relation holds 


Upp = ul (11) 


Setting in this equation f =2 s and extracting the root of the 


n-th power, we obtain 


Next setting in equation (11) t= 1 s and substituting p for n, we 
obtain 


u, = ut 
It follows from the last two equations that 
2 
1 = ut 
n 


Thus, for any positive rational number ¢ (i.e. the number of 


the form 4, where p and n are positive integers) it holds 


u, = ui, 
or returning to the previous notation 


E(u 
ae ee cee (12) 


Here v, is the velocity of the falling body at the instant t= 1 s. 
Inasmuch as relation (12) is valid for rational t’s it is valid for 
all t's. 
Let us take, for example, the instant t=//2 s = 1.414... s. 
The numbers 1.4, 1.41, 1.414, etc. are rational and for this reason 
relation (12) holds for these values of t: 


1.4 
M4 vy — 
oe \a® 
1.41 
Daa — ler ee 
= Sate (13 
», — 28. >, — 2. 
0 b 0 b 


Should we set ¢ equal to rational values each of which is 
progressively a more exact approximation of the number y2 
the left-hand side of equation (13) would approach the limit 


y2 
Vie ae 
——_——— and the right-hand side the limit ” 
aaa F oe. 
‘0 b Vo b 


* See, for example, the concept of an irrational power in the course 
of algebra by A. P. Kiselev, part II (p. 97 in the 1960 Russian edition). 
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Thus we obtain in the limit 


mg mg \V? 
Yi b ar 
i 
Oo % OE 


Similar treatment is, of course, applicable not only to va but 
to any irrational value of t. Hence, relation (12) is valid for any t. 
Introducing the notation 


=c 
we obtain from (12) 
», — md 
b _¢ 
mg 
ae 
whence 
o= + (0-H )e (14) 


Formula (14) obtained above for the velocity of a falling body 
is not the final one since we do not know the value of the number 
c which enters this formula. To calculate the number c, let us 
find from formula (14) the acceleration of the body at the initial 
instant of its motion. The average acceleration for the first h seconds 
has according to (14) the following value 


ig a 
4 _— Una %o _ B+ (v6 2) ae _ mg ch 1 
av h =, h =| Uo b h 


As h tends to zero, this expression gives us the acceleration dp 
at the initial instant: 


‘ F ch 1 
to= Jima = fm(oo— EEE a 
Denoting the expression c’ — 1 by x, we obtain 
ct 1 =x, ch=l4+-x, h = log, (1 + x) 
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mg \ c'—1 


Thus instead of the expression { v9 — —> which enters 


b h 
(15) under the limit sign, we obtain the expression 
co et IND: _ ng 
mg x oe jas 
’o~ “h} Tog. (1+) 


1 7 ae 
= loge + x) log. (1 + x)* 


Note that as h tends to zero the limit of the number c" is 
unity and that the number x = c" — 1 approaches zero. Hence we 
may write: 


mg 
“ae Obs. 
ao = lim — ere (16) 
log, (1 + x)* 


1 

The limit of the expression (1 + x)* as x tends to zero is termed 

the number e. We shall leave without proof that this limit exists, 
1 

i.e. that the expression (1 + x)* does in fact approach some value 
as x > 0. The proof (by the way, quite elementary) of the existence 
of this limit may be found in the first chapters of any textbook 
on higher mathematics. * 1 

We shall confine ourselves to calculating the values (1 + x)* for 
x=0.1, x=001, x =0.001, and x = 0.0001. The results of such 
calculations are presented below (these calculations may be 
performed with the aid of logarithmic tables, preferably seven- 
place ones; one may also use Newton’s binomial formula): 


ate 


(14+0.1)°° = 1.19 ~ 2.59374 
(1+ 0.0107 = 1,011 ~ 2.70481 
(1+ 0.001) 007 = 1.00119 ~ 2.71692 
(1+ 0.0001) 20001 = 1,0001!00 ~ 2.71814 


*) See, for instance, V I. Smirnov, Kurs vysshei matematiki (A Course 
in Higher Mathematics), V. 1 
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These calculations clearly demonstrate that as x0 the limit 
1 


of the expression (1 + x)* is e = 2.71... . We now obtain from (16) 
5, — 2. 
0 


es 
log, e 


On the other hand, it follows from (5) that 


b 
ao =~ (%--) 


a = 


Tb bmg 
loge j.m\° 5 
whence 
_m _> 
loge = — e=c 4, c=e ™ 


-+ 
Finally substituting e ™ for c into formula (14) we obtain formula 
(7) thus concluding the proof. 


Differentiation 


The Concept of the Derivative 


Equation (5) as we have just seen may be solved quite accurately. 
This equation relates the quantity v (the velocity of a falling body) 
with the quantity a which shows how rapidly the quantity v changes 
(acceleration is “the rate of change of speed”). 

When we talk of the rate of change of some quantity, we assume 
that the quantity we are dealing with is not a constant characterized 
by a number, but a variable quantity, i.e. a quantity whose value 
changes with time. Here are some examples of quantities that 
change with time (or depend on time): the velocity and the accelera- 
tion of non-uniform motion, the magnitude of alternating current, 
ete. 

Let y be a quantity whose value changes with time. For the sake 
of convenience we shall denote by y, the value this quantity 
assumes ¢ seconds after the process under observation has started. 
The difference y,., — y, shows how much the variable y has changed 
in h seconds (in the interval t, t + h seconds after the process has 
started). The ratio 


ae Ne (17) 
for its part shows how much y changes on the average every 
second (within this interval of time), i.e. this ratio represents 
the average rate of change of the variable y. Putting h progressively 
smaller, we shall obtain the values of the average rate over 
progressively shorter intervals starting from the instant t. As a limit- 
ing case (when h tends to zero) relation (17) yields the rate of 
change of the variable y at the instant t. We know already that 
this rate of change is designated in mathematics by 


lik V+h— Ve (18) 
h-0 h 
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Expression (18) is termed the derivative of the quantity y with 
respect to time t, and as we have seen it represents the rate 
of change of the variable y. We may deal with a variable which 
changes not with time, but depends on some other quantity. 
For instance, the area of a circle depends on its radius. Denoting 
the area of a circle of radius R by Sp, we obtain 


Sr = mR? (19) 
Should we investigate the dependence of the circle’s area on its 

. ‘ A Sr +h Sr ; 
radius we would arrive at the ratio aaa ae which expresses 


the mean rate of change of the area with the change of radius. 
The limit of this ratio (as h-0) is the derivative of S with 
respect to R. 

The concept of the derivative is one of the basic concepts of 
higher mathematics. In cases when the variable y changes in 
accordance with variations of a quantity x (or when, as may also 
be said, y is a function of x) the derivative of y with respect 
to x is denoted by y’, or more often, by one of the symbols 
dy ad 


dx’ dx?’ 


dy = lim 22th x 


here the letter d should not be cancelled out since it is not 
a multiplicand but denotes the process of finding the derivative, 
or, as is also said, the differentiation. 


Let us by way of example calculate the derivative & of func- 
tion (19): 
dS lim Sreh = Sr li 1 (R + h)? = mR? 
— = 1 TO - ——————— A ———— I, “ort, 
dR woh ho h 
= lim (2nR + th) = 2nR 
-0 
i.g. the derivative of the circle’s area with respect to its radius 
equals its circumference. 
For the second example let us find the derivative ee of the 


dt 
distance with respect to time. Denote the distance covered by 
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a body up to the instant t (i. e. t seconds after it has started moving) 
St+n — St 
h 
during the interval t, t + h, the limit of this ratio for h > 0 being 

the value of the speed at the instant t: 


by s,. Then the ratio will represent the average speed 


li Stth — St ds 
v, = lim ————— = —— 
: h-0 h dt 
. .. _ dv 
In the same way one may calculate the derivative Te The ra- 
tio 
Drth — Uf 
h 


is the average acceleration during the interval t, t+ h, the limit 
of this ratio being the value of the acceleration at the instant rf 
(compare with what was said on p. 16): 


ys, Orth Ur dv 
gt hat 
The relations that have just been proved 
ds 
v= "Th (20) 
and 
dv 


play an extremely important part in mechanics. 


The Differential Equation 


Let us again turn to equation (5). According to (21) this equation 
may be rewritten in the form 


dv b mg 
=— - 22 
dt ~~ m\°— 6 oa 
Now it is evident that this is an equation with one unknown 
quantity v. However, this is not an algebraic equation, but an 
equation relating the quantity v to its derivative. Such an equation 
is termed a differential equation. Comparing differential equation 
mg 
b 


(22) with its solution (7) and denoting a by k and by c, we 


may make the following statement. 
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THEOREM. The solution of the differential equation 


dv 
a —k(v —c) (23) 
is represented by the expression 
v=c+(v9—cde™ (24) 


where vo is the initial value of v (i.e. its value for t = 0). 
In future we shall be able, using this theorem, to calculate some 
other physical phenomena, as well. 


Two Problems 
Leading to Differential Equations 


(a) Switching on an Electric Current. Consider an_ electric 
circuit consisting of a coil and a battery (Fig. 2). The electrical 


Fig. 2 


properties of a coil are rather complex, but for some cases they 
may be described with great accuracy with the aid of two quantities: 
its resistance and its inductance. The coil is in fact usually represented 
by two sections connected in series: by the resistance and the 
inductance (Fig. 3). The voltage drop across the resistance is 
proportional to the current i passing through the coil (Ohm’s Law) 


V=Ri 
28 


Coil 


L 


Switch 


Ws a 


The proportionality factor R is termed the resistance of the coil. 
The voltage drop across the inductance is proportional to the rate 
of current variation. Denoting the current variation rate by w 
(it is measured, for example, in amperes per second) and the propor- 
tionality factor by L, we obtain for the voltage drop the expression 


V =Lw 


The quantity Lis termed the inductance of the coil. The voltage 
drop across the coil is the sum of voltage drops across the 
resistance and the inductance, i.e. it is expressed by the formula 


V=Lw+ Ri (25) 
Formula (25) is well supported by experiment (if the frequency 
of the current passing through the coil is not too high). This is 
the formula we are going to use. Denote the electromotive force 
(emf) of the battery by E. Equating the emf of the battery to the 
voltage drop across the coil, we obtain on the basis of Kirchhoff’s 
Second Law the equation (we neglect the internal resistance of the 
battery and the resistance of the conductors): 


E=Lw + Ri 


w= -F(i-$) (26) 


The solution of this equation may be easily obtained from 
the theorem formulated on p. 28. Indeed, denoting the current at 
the instant ¢ by i,, we may contend that the quantity 


or 


ban —h 


Way h 


is the average rate of change in the current during the interval 
t, t+h. When h-0 we obtain the rate of the current change 
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at the instant t 
daa — i di 
w = lim beth 7 he ae 
h—0 h dt 
Hence the quantity w is the derivative of the current i, and 
equation (26) may be rewritten in the form 


The only difference between this equation and equation (23) is 
that the variable we are seeking is denoted by i and not by v 
which is, of course, inessential. The constants k and c contained 
in equation (23) in this case assume the following values 

R E 
Necgeh Rg 
Hence the solution of this differential equation will be of the form 


(see (24)): 
he i ge a 
Por oT R 


Should the current ip the instant the battery was switched on 
(t = 0) be zero the formula would assume a simpler form 


i=l opt 


It follows from this formula that the current, zero at the 
instant of switching on, grows continuously approaching the value 


* i.e. the value of the current which would pass through the 


coil if it had the same resistance R and no inductance at all. 

(b) Radioactive Decay. Suppose we have a piece of rock 
which contains some radioactive material. The atoms of the radio- 
active material decay and are transformed into another chemical 
substance — the product of decay. Hence, in the course of time 
the amount of the radioactive material contained in the piece 
of rock decreases. Let us introduce the concept of the decay rate. 
Suppose at some instant t the amount of the radioactive material 
contained in the rock had been m, grammes and h years later 
it decreased (because of the decay) to m,,, g. The expression 


M+h— ™ 
h 
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shows by what amount, in grammes per year, the mass of the 
radioactive material has on the average decreased (during the 
respective period of time). It would be natural to term this 
expression the average (mean) decay rate for the period. The limit 
to which this expression tends as h- 0 is the decay rate at the 
instant t. Denote it by u 


aS iim ee 2 ae 
~~ hao h dt 


Note that the decay rate is negative since the mass of radioactive 
material decreases with time. 

What does the decay rate depend on? When the amount of 
radioactive material contained in the rock is small, it may be 
accepted that the decay rate is proportional to the amount of radio- 
active material contained in the piece at a certain instant, i.e. 


u=—km 


where m is the mass of radioactive material and k a positive 
constant (proportionality factor). 

The arguments in favour of approximate validity of this law are 
simple if one assumes that the decay of some atoms of radioactive 
material does not affect the state of the others. In these conditions 
it may be assumed that the amount of radioactive material decaying 
per unit time out of every gramme is always approximately the 
same, say k grammes, no matter how much radioactive material 
still remains in the rock. In this case the amount of radioactive 
material decaying per unit time out of m g will be km g. 

Is the assumption that the decay does not affect the state of 
the remaining radioactive atoms correct? One should bear in mind 
that the decay products may strike another atom, thereby causing 
it to decay; new particles will thus be generated — this may lead 
to the decay of further atoms, etc. Such a chain reaction (this, 
for example, is the basic process in the functioning of the atom 
bomb) would not fit with the assumption of the independent decay 
of the atoms. To make this chain impossible, the particles emitted 
in the course of the decay should be lost and should not (in most 
cases) strike other radioactive atoms. This will be the case if 
radioactive material makes up only a small percentage of the mass, 
the bulk of the rock being non-radioactive. Then the overwhelming 
majority of the particles emitted in the course of the decay will 
be lost in collisions with non-radioactive atoms of the rock, and 
a chain reaction will not be possible. Therefore, for small amounts 
of radioactive material contained in the rock the decay of the 
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individual atoms may, to a good approximation, be assumed to be 
independent. 

We obtain accordingly the following differential equation for 
evaluating the mass of the remaining radioactive material 


which differs from equation (23) in that the unknown quantity 
is designated by the letter m instead of v, and that the constant 
c is in this particular case zero. Hence, according to (24) the 
solution will take the form 
m, = moe" (27) 

where my is the mass of the radioactive material at the initial 
instant (when we begin to take interest in the process of decay). 

Example 5. How many years will it take for the amount of 
the radioactive material to decrease by half? 

Solution. To obtain an answer to this problem, we should 
determine ¢ from the equation 
-kt 


Moe "= > Mo 
Cancelling out mo and taking the logarithm, we find 
1 0.69 


ie ace 


This period is termed half-life of the corresponding radioactive 
material. Note that this period does not depend on the initial 
amount of the radioactive material, but only on k, i.e. on the 
nature of the material. For instance, half-life of radium is 
1590 years, that of uranium 238 — about 4.5 thousand million years, 

Example 6. Formula (27) permits us to draw certain conclusions 
concerning the age of the Earth. 

Suppose a piece of rock extracted from the depths of the Earth 
contains together with some inclusions m g of the radioactive 
material and p g of its decay product. Suppose in addition that 
every gramme of this radioactive material yields (upon complete 
decay) r g of the decay product. This means that p g of the decay 


product was produced from 4 g of the radioactive material. 


Hence if we assume that the radioactive decay process started at 
some instant (i. e. up to that instant the material had not contained 
a single atom of the decay product), we must put the initial mass 


of the radioactive material equal to m+ £ To calculate the time 


that has passed from that imaginary instant (since the beginning 
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of the decay) until today, we must in accordance with (27) solve 


the equation 
m= (m + 2) e7* 
_ r 


with respect to t, whence we obtain 


1 P 
t= {log (1 + £) 


Such calculations performed for some minerals contained in the 
Earth yield approximately the same value for tf — around 
2x 10° years. Thus the conditions on the Earth that permitted 
a normal decay process persisted for several thousand million years. 
It may be conjectured that some thousand million years ago 
the matter which now makes up the Earth existed under quite 
different conditions which facilitated the production of radioactive 
atoms from the simpler atoms and from other particles. 


Napierian Logarithms 


The formulas for solving the problems presented above contain 
the exponential function with the base e. When making calcula- 
tions based on these formulas with the aid of logarithmic tables, 
one can exclude some operations if one uses logarithms with the 
base e. For instance, taking the logarithm in formula (27) with 
the base e and with the base 10, we obtain 


log, m, = —kt +log,mo, logio m, = —kt logio e + logio mo 


In the second instance an extra logarithm has to be taken and 
an extra multiplication performed. Moreover, problems result 
in forrnulas which contain logarithms with the base e as we have 
seen in Examples 5 and 6. The number e comes up in other 
mathematical problems as well, and it is very convenient to use 
logarithms with the base e, espccially in theoretical questions. 
Logarithms with the base e are termed Napierian, or natural, loga- 
rithms and denoted by the symbol In: the meaning of the expression 
In x is the same as that of log, x. Common and natural logarithms 
are related by the expression 


logio xX = M-Inx 


where 
M = logio e ® 0.4343 


This relation is easily obtained if one takes the logarithm with 
the base 10 of the identity e"* =x 
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Harmonic Oscillations 


The Problem of Smal} Oscillations 
of a Pendulum 


Suppose a thread of length | is fixed at point C with a body 
M suspended from its other end (a pendulum). The problem is 
to find out the nature of the motion of the body M. To solve 
this problem mathematically, we shall simplify it somewhat. First 
of all we shall assume that the thread from which the pendulum 
M is suspended is inextensible and weightless. 

We shall examine the motion of the pendulum which takes place 
in a single vertical plane passing through the suspension point. 
Because the thread is inextensible we may be sure that the body 
M will move along a circumference of radius | with the centre 
at point C. The assumption about the negligible mass of the thread 
means that its weight is negligible as compared with that of the 
body M. This makes it possible for us to assume that external 
forces act only on the body M. We shall assume the pendulum 
M to be a heavy point (i.e. we shall attribute to it a mass m 
and neglect its dimensions). From all the forces that act on the 
body M we shall in addition to the tension of the thread take 
into account only the force of gravity. When considering this 
problem, we may also neglect the air resistance (for instance, we 
may presume the body M to be confined to a hermetic evacuated 
vessel). The problem of how much the motion of the pendulum 
in the atmosphere differs from its motion in a vacuum is treated 
in the note on p. 49. 

Let the body M at a certain instant be at point A on the 
circumference along which it is moving. Denote the lowermost 
point of the circumference by Q, the length of arc QA by s, and 
the central angle 2 QCA corresponding to it (Fig. 4) by a 
(in radians). Then 


s=lo (28) 


We shall assume the arc s and the angle a to be positive when 
point A is to the right of Q and negative in the opposite case. 
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Fig 4 


Now let us deduce the equation from which we are to learn 
the law of the pendulum’s motion. The difference in altitude of 
points A and Q is 


h=QB = CQ—CB=I1-—Icosa=I(1 — cos a) = I 2 sin? 2 
Assuming the potential energy of the pendulum to be zero at 

point Q, we find for the potential energy of the pendulum at point 

A the value 

Bae 


W”) = mgh = mg2l sin? 5 


The kinetic energy of the pendulum is 


we — mv? 
2 
where v is the velocity of the body M. 
Hence the total energy E of the pendulum (at point 4A) is 
expressed by the formula 


2 
E =e + Qmgl sin? (29) 
2 2 

The pendulum does not perform any work in its motion (we 
neglected friction and drag!), and for this reason its energy E 
remains all the time constant. 

We shall somewhat simplify equation (29) by treating the 
problem of small oscillations of the pendulum, i.e. of such motion 
of the pendulum in the course of which the deflection angle from 
its equilibrium position Q remains small. Let us clarify what we 
mean by “small angles”. The point is that an exact solution of 
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equation (29) cannot be written with the aid of any known opera- 
tions. Therefore the question arises as to whether equation (29) 
may be replaced by a simpler one. Naturally this simplification 
should be such that the solution of the simplified equation would 
be a very accurate approximation to the solution of equation (29). 
Note that hereby we do not in principle introduce a new inaccuracy 
since relation (29) is in itself approximate*’, so the applicability 
of an approximation of a definite sort depends only on the degree 
of accuracy we want to obtain. 

The usual simplification to which equation (29) is subjected is 
the substitution of an angle @ for sin @. That such a substitution 
is permissible for small angles @ will be seen from Fig. 5, which 


Fig 5 


shows an arc A’Q’B’ of the circumference of radius C’Q’ = 1; 
the angle @ is plotted on both sides of radius C’Q’. The length 
of section A’B’ is 2sin@ (for A’S’ is the sine line), the length 
of arc A’B’ being 29 (the angles are, of course, measured in radians!). 
It is clear from the figure that small angles @ differ little, this 
difference being the less the smaller is . For instance, it may 
readily be checked with the aid of trigonometric tables that 
for angles not exceeding 0.245 radian (i.e. ~ 14°) the ratio 


* In deducing equation (29) we have made some simplifications: 
neglected the drag, the weight of the thread, the dimensions of the 
body M, etc. 

Note that any physical law or any mathematical relation linking physical 
quantities (for example, relations (1), (2), (3), (4), (5), (25), and (29)) is 
an approximate one, for in every case there are “forces” which were not 
taken into account when the physical law or mathematical relation was 
formulated. This, however, does in no way belittle the enormous importance 
of physical laws. For instance, the accuracy of Newton’s Second Law 
and of Ohm’s Law is under ordinary conditions extremely great. 
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sin p 


differs from unity by less than 0.01; for angles less than 

1° (0.017 radian) this difference drops to about 0.0005. 
Accordingly, assuming the deflection angle of the pendulum 

to be small, we substitute = for sin = i.e. we substitute a new 


equation for equation (29) which differs little from the former: 
mv? a \? 
> + 2mgl (5) =E 
Taking into account (28), we may write this equation in the form 
mv? mgs? 
2 +3 


ae eee 
g 


=E 


or 


= (30) 


This equation contains two unknown variables: s and v (we 
assume that constants g, |, m, and E are known). Yet this equation 
(same as equation (5)) may be solved, for quantities s and v are 
not arbitrary but related by expression (20). It follows from (20) 
that equation (30) may be written in the form 


I (ds\? |, 2E 
5(%) +S = mg. (31) 
so it is, in effect, an equation of one unknown. Let us solve this 


equation. 
Let us choose a coordinate system lying in a plane with the 


value of s measured along the x--and the value of \t v along 


the y-axes. At any instant ¢ there is a definite value of the distance 
s and the velocity v of a body M, i.e. a definite point N in the 
plane (Fig. 6). Conversely, if we know where point N lies, we can 


Fig 6 
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; . i ; i 
find its coordinates s and \t v, i.e. we can learn the position 


of the pendulum and its velocity. Thus at any instant ¢ the 
pendulum M is symbolized by some point N. The length of section 
ON is easily calculated with the aid of the Pythagorean theorem: 


ON =|/PN? + OP? = lie +3? 


i.e. (in accordance with relation (30)) 


on = |/2# 

mg 
As the pendulum swings, the values s and v will change, i.e. 
point N will move across the plane which contains the coordinate 
system. But the distance from point N to the origin of coordinates 


F , 2lE : 
always remains the same, i.e. equal to |/——. Hence point N 


moves along a circumference of the radius 

2IE 

mg 

This circumference is termed phase circumference. 
Let us find the velocity with which point N moves along 

the circumference. This velocity is directed along a tangent to the 

circumference. Let it be represented, for example, by a vector NA 

(Fig. 7). Let us decompose the vector NA into a horizontal and 

a vertical components. Then a horizontal component NB will 

represent the velocity of motion of point P along the x-axis. 

The distance from point P to O being s, the velocity of motion 


R= (32) 


of point P is Fao, ie. NB =v. 
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Now since triangles ONP and NAB are similar, we obtain 
PN:ON=NB:NA_ or besa ona 
From the latter proportion we obtain 


NA=R|/4 
Such is the velocity with which point N moves along the circum- 
ference. 
Denote the initial deflection and the initial velocity of a pendulum 
by so and vo, respectively, and the corresponding point on the 
phase circumference by No. Then the radius of the phase circum- 


ference will be 
R= \z ve + sz (33) 


(see (30) and (32)). An angle gp = 2 XON, will be determined 
from the relation 

sy 
g 0 
So 


tan Op = (34) 


(Fig. 8). Next, t seconds after the motion of the pendulum has 


Fig 8 
started, point N moving at a velocity R \/s will travel a distance 
NoN=R se: along the phase circumference, and an angle 
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Zz NoON will accordingly be ye As a result (Fig. 9) 


Fig. 9 


0= £XON = 2 XONy~ <.NON = oo — |/¥ 


Hence we obtain 


OP = Roos @ = Reos (9 — 41) =Reos(|/21- 0) 
— Remm — Re g9\_ : ] 
PN = Rsing = Rein ( - 41) = -Rsin(|/41— 0) 


Remembering finally that OP = s and PN = yi0 we obtain 


s= Reos(|/41—e). --|/4 Rsin(|/4- 00) (35) 


These formulas express the deflection and the velocity of the 
pendulum ¢ seconds after it has started to swing, i. e. they represent 
a complete solution of problem of pendulum’s motion (with the 
simplifications we have made). Let us consider some examples. 

Example 7. The pendulum is initially deflected to the right 
to a distance sy and dropped with a zero initial velocity. Find 
its deflection and velocity at the instant t. 

Solution. In this case R = So, ®) =0, and from formulas (35) 


we obtain 
s= sycos|/& on -|/S usin |/ 44 


Example 8. The pendulum was initially in the equilibrium 
position Q and was thrown out of this position by a thrust which 
imparted to it the initial velocity vg directed to the right (i.e. 
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a positive initial velocity). Find the deflection and the velocity 
at the instant ¢. 
Solution. From formulas (33) and (34) we obtain in this case 


that R= Jie o =>: and from formulas (35): 


s= [iam(/-8 = |Eeosin|/4 
9, 


b= —Uo0 sin i -)- Uo COS ee 


Example 9. Let us find the derivatives of functions sin wt and 
COS wt. 

Since v is the time derivative of s, comparing the values of s 
and v in Example 8, we conclude that 


d I : g.\_ g 
A (Yqeeosin |/%1) = vcos |/% 


Similarly, from Example 7 we obtain 


$e (wom |/$0) =~ /$ sonia [4 


Setting, for instance, vp = \s So = 1 and denoting the quantity 


\3 by , we obtain from these formulas: 


d ., d ; 
rr Sin Wt = @ COS wt, ih cos ot = — sin wt (36) 


Example 10. Since cosine and sine are periodic functions 
after some time T, termed the period of oscillations, the pendulum 
again returns to the initial position and again repeats the same 
motion. Let us find the period of oscillations of the pendulum. 

Solution. The values of sine and cosine remain unchanged when 
their arguments change by 2x. Therefore the period of oscillations 
of the pendulum is the period T during which the expression 
in parentheses following the sine and cosine symbols in equations 
(35) increases by 2n. In other words, the difference between the 


values of the expression \s: — Qo at the instants t and t+ T 
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must be 22 


|/$0+n-o=(|/4r—o0) +2 


From this relation T is easily obtained 


T =2n yt (37) 
g 


In this way the motion is periodically repeated every T seconds. 
The pendulum oscillates periodically. During each period (i. e. every 
T seconds) the pendulum, as may be seen from (35), once assumes 
the extreme right-hand position (cosine becomes equal to +1) 
and once the extreme left-hand position (cosine becomes equal 
to —1). At these instants of maximum deflection the velocity 
of the pendulum turns zero (see (35)). When cosine turns +1, 
the sine of the same argument turns zero. The maximum velocity 
of the pendulum is attained (sine turns +1) when it passes through 
point Q (cosine turns zero). 


The Differential Equation 
of Harmonic Oscillations 


We have deduced formulas describing the motion of the pendulum 
from equation (30) or, which is the same, differential equation (31). 
There is another differential equation which, too, describes the 
motion of the pendulum we have been discussing above. Its deduc- 
tion is very simple. 

Suppose a body M at some instant is at point A of the 
circumference along which it moves. Let us decompose the force 
of gravity (which we assume to be equal to mg and to be directed 
vertically downwards), using the parallelogram rule, into two 
components: one tangent to the circumference at point A and 
the other perpendicular to this tangent. The component which is 
perpendicular to the tangent tends to stretch the thread and is 
compensated by the tension of the thread (for the thread is presumed 
to be inextensible). A force F acts along the tangent in the direc- 
tion of point Q and is equal in value, as may readily be seen, 
to mg sina (Fig. 10), i.e. the force F is negative for positive 
a’s, and vice versa. Thus 


F = —mg sina 


Tf both mutually compensating forces — the tension of the thread 
and the component of the force of gravity perpendicular to the 
tangent — are ignored, there remains only the force F acting on the 
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Fig 10 


body M (we neglected the air resistance), and we may therefore 
write on the basis of Newton’s Second Law 


ma = —mg sin a 
or 
a= —gsina 
Recalling that we are now dealing only with small oscillations 
of the pendulum, and that this makes it possible to approximate 
sina by a, we may write this equation in the form 


= —ge 


or according to (28) 
at 4 s=0 (38) 


This is just the required equation. Let us show how it may be written 
in the form of a differential equation. It follows from relations 


_ dv 


ds ; ater ; 
a= and v =—— that if we take a derivative of a distance s 


dt dt 

once and then take the second derivative of the first one obtained 
(i.e. of the velocity), we shall obtain the acceleration. In other 
words, the acceleration is the second derivative of the distance s 
(with respect to time t). This is written down as follows 


_4 (as 
a= at \ at 


a=—> (39) 


or 
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2, 
The symbol “s (the second derivative of s with respect to 1) is 


not treated as an algebraic expression, but as a single symbol; 
we are not allowed to perform any operations with its components 
(for instance, to cancel anything out of the “fraction”). It follows 
from (39) that equation (38) may be written in the form of 
a differential equation as follows 


d? 
Sr +is=0 (40) 


Note that we are already able to solve this equation. Indeed, 
it describes the law of variation of the quantity s, i.e. the law 
of the oscillations of the pendulum, and this we have already 
investigated. Therefore we may state at once that the solution 
of equation (40) is contained in the first formula of (35). We may 
express this in more detail as follows. 

The solution for differential equation (40) is 


s= Roos (|/41 - 0.) 


where R and @o are given by formulas (33) and (34). Note that 
to find the quantities R and @) we must know the initial deflection 
So and the initial velocity vp, i.e. the values of s and 
& taken at the initial instant. 


Denote a quantity \s by . Then we may express the state- 


ment made above in the following form. 

THEOREM. The differential equation 
es + @’s =0 41 
dt? _— ( ) 


has the solution 


where R and @o depend on the initial values of s and & 
t 

Equation (41) is known as the equation of harmonic oscillations. 
Any quantity described by this equation is said to be oscillating 
harmonically; this means that the respective quantity changes with 
time in accordance with the law (42). The quantity @ which enters 
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differential equation (41) and its solution (42) is termed angular 
oscillation frequency, and T = = oscillation period. When the 


quantity s oscillates harmonically, any value it assumes is repeated 
every T seconds (see Example 10). 

Let us compare differential equations (23) and (41). Equation (23) 
contains only the first derivative, therefore it is termed a first-order 
differential equation. Equation (41) is a second-order differential 
equation because it contains the second derivative. Note that only 
the initial value of the quantity v itself had to be known to solve 
the first-order equation (23). To solve the second-order equation 
(41), on the other hand, we must know not only the initial value 


of the quantity s itself, but of its derivative & as well. In short, 


to solve a first-order equation we must know the initial value of 
only one quantity; to solve a second-order equation, of two. 

Note that we have obtained the solution of equation (40) from 
physical considerations: both equations (31) and (40) describe the 
same physical phenomenon, and for this reason they should have 
identical solutions expressing the law of pendulum’s oscillations. 
Of course, such argumentation is guesswork and not a rigorous 
mathematical proof. It may be proved by purely mathematical 
means that equations (31) and (40) are equivalent, i.e. that their 
solutions are the same: after differentiating both sides of equation 
(31) we obtain equation (40). Conversely, one may obtain equation 
(31) from equation (40), but then one will have to use an operation 
reverse to that of differentiation. Such an operation (termed 
integration) forms together with differentiation the basis of higher 
mathematics as a whole. A more detailed discussion is beyond 
the scope of this booklet. 

Still, by applying formulas (36), the reader will easily find the 
second derivative of function (42) and make sure for himself that 
the function satisfies equation (41). 

Let us consider two examples from the field of physics which 
lead to the equation of harmonic oscillations. 


The Oscillatory Circuit 
Let us consider an oscillatory circuit, i. e. a closed electric circuit 
consisting of a coil and a capacitor. The coil possesses an inductance 
and a resistance (see p. 29). The complete circuit may be represented 
in the form of a circuit diagram (Fig. 11). Let us denote by q 
the electric charge which goes over from one plate of the capacitor 
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Coil 


t 


Switch 


Capacitar 


Fig 11 


to another and by i the current in the circuit. (We presume that 
initially the capacitor had a charge qo, and that a current ip was 
in the coil, and are interested in subsequent changes in these 
quantities.) In that case the voltage across the capacitor will be 
4 where C is its capacitance; the voltage drop across the coil 
will be equal to Lw+ Ri, where R is the resistance and L the 
inductance (see (25)). According to Kirchhoffs Second Law the sum 
of the voltage drops in the circuit is zero, i.e. 


Lw+Ri+4=0 (43) 


A quantity i is the derivative of q with respect to t. Indeed, 
if the values of the charge q were q, and q,,, at the instants 
t and t+h, respectively, this means that during that interval 
a charge passed through the cross section of the circuit’s wire 
(anywhere in the circuit) equal to q,.,— 4, Therefore the average 
current over the interval t, t+ h was 


pias Q+h— 
av h 
Taking the limit, we obtain 
. - Gt4+n— & dq 
= lim ———— = — 
: 436 h dt 
: . dq di 
It follows from relations i = ap and w= Wi that w is the de- 
rivative of i= a i.e. that w is the second derivative of q: 
dq 
wane <a 
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Hence equation (43) may be rewritten in the form: 
d 
Ai Roy 20 (44) 
The differential equation we have just obtained is more complex 
than equation (41), since in addition to the unknown function q 


: Sankar ge 3 : etd 
and its second derivative <r it contains also the first derivative 


dq 
dt 
note on p. 49), but shall only discuss the case when the resistance 
of the coil R is negligible (compared to the quantities L and C), 


. We shall not, however, try to solve equation (44) (see the 


: ; dq . : 
neglecting on this account the term R < in equation (44). Then 


the equation will assume the form: 


dq 4 
tae a 
or 
2 


Equation (45) is, evidently, an equation of harmonic oscillations 
(see (41)) with the oscillation frequency of the circuit 


1 
QO= —— 
VLC 
and the period of oscillations expressed by the formula 
T=2n VLc 


The solution of equation (45) is of the form (see (42)): 


4= Roos (Hi 0) 


where R and @po depend on the initial values, i.e. on go and ipo. 


Oscillations Produced by 
the Action of the Elastic Force 
of a Spring 
Suppose a weight of mass m is suspended from a spring. 
The spring will be somewhat elongated under the action of the 
force of gravity (the elongation will stop as soon as the spring 
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tension compensates the force of gravity); in this position the weight 
and the spring can remain stationary (be in a state of equilibrium). 
Should we shift the weight from its equilibrium position pulling 
it downwards, the spring tension would exceed the force of gravity 
and the resultant force would be directed upwards. Should we, 
on the contrary, shift the weight to a point above the equilibrium 
position, the resultant force would be directed downwards. Thus, 
this resultant force “tries” to return the weight to its equilibrium 
position. 

For the sake of simplicity we shall content ourselves with only 
the vertical motion of the weight: up and down. Denote the posi- 
tion of equilibrium by O, the position of the weight at some 
instant by A, and a distance OA by s, assuming the positive 
direction along the vertical to be downward from point 0, i.e. 
we will take s to be positive if the weight (point A) is below point 
O, and negative if it is above point O. Denote the resultant force 
of the spring’s gravity and tension by F, and the air resistance 
by S. We shall assume that no other forces except F and S act 
on the weight. In accordance with Newton’s Second Law we may 
write: 


ma=F+S 


where a is the weight’s acceleration. The force F which tends 
to return the weight to its equilibrium position is directly propor- 
tional to a deviation s, i.e. it is numerically equal to the quantity 
ks, where k is a proportionality factor. This assumption is well 
supported by experiment (for not very great deviations from the 
equilibrium position). The quantity k is termed the stiffness of the 
spring. For positive s (point A is below O) the force F is directed 
upwards, i.e. it is negative, for negative s the force F is positive. 
In other words, the sign of the force F is opposite to that of the 
deviation s, i.e. 


= —ks 


We shall assume for the force S the same value as above (see 
(3)), ie. 


S= —bv 


In this way we obtain the following equation of motion of the 
weight 


ma = —ks — bv 
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or 


ma+ bv +ks=0 (46) 
Since v = é and a= Lee this equation may be written in the 
form 

mos 4 ods +ks =0 (47) 


Differential equation (47) is analogous to differential equation 
(44), which we obtained while solving the problem of the oscillatory 
circuit. We shall not solve equation (47) (see the note below), 
but consider only the case when the air resistance may be neglected 
(i. e. when a quantity b is very small compared to quantities m 
and k). In that case equation (47) takes the form 

ds ek 


Relation (48) is an equation of harmonic oscillations with the 


frequency = 
k 
VE 


= 2n \z 


The solution of equation (48) is according to (42) of the form: 


/* 
s= Reos( 5 t-¢) (49) 


where R and @p depend on the initial conditions, i.e. on so and v9. 

Note. In order to obtain the equation of harmonic oscillations 
we neglected the forces of friction and air resistance when studying 
the oscillations of a pendulum and a weight suspended from the 
spring, and the resistance of the coil when studying the oscillatory 
circuit. Physically it means that no energy is spent in the process 
of the oscillations, mathematically this leads to the omission of 
the term containing the first derivative. As a result we have obtained 
harmonic oscillations, i. e. oscillations that remain constant from 
cycle to cycle, or sustained oscillations. 

What difference would it make if while solving the above 
problems we took into account the force of air resistance and 


and the period 
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the voltage drop across the resistance? For instance, what is the 
difference between the solutions of equation (44) and those of equa- 
tion (45)? Mathematical calculation (which we shall not present 
here) shows that when R is not too great, equation (44), too, 
describes an oscillatory process. The oscillations described by equa- 
tion (44), however, become weaker with time; for this reason they 
are termed damped oscillations. The physical explanation of this 
difference is that the energy of oscillations constantly decreases 
being transformed into heat because current passing through the 
resistance R generates heat. In the same way the oscillations 
of the pendulum, too, become gradually weaker because due to 
friction and air resistance the pendulum’s energy is gradually spent 
on heating the pendulum itself and the surrounding atmosphere. 
However, when the resistance is small, the difference between the 
damped and the sustained (harmonic) oscillations during a short 
interval of time (for example, for several periods) is not great. 
The damping makes itself felt only after a substantial time. Should 
we, for example, suspend a heavy weight from a string and deflect 
it slightly from its equilibrium position, the decrease in the amplitude 
of the oscillations after 10-15 periods would be quite negligible 
and unnoticeable to the eye. We would be able to observe it only 
several minutes after the oscillations have started. 

Here for the purpose of comparison is the exact solution of 
equation (47) (without deduction). We shall presume that the value 
of the factor b in the expression_for the force of air resistance 
is not too great (namely, b < 2|/ mk). Then the solution of equation 
(47) takes the form 


b, Tb 
s=Re ™ cos (/* - () t— a) (50) 


where R and @o are determined by the initial conditions. It may 
be seen from this formula that s decreases infinitely with time 


_ 
(multiplicand e 2m‘ becomes less and less as time ¢ increases). 
Figure 12a, b shows the graphs of functions (50) for various values 
of a factor ~—- The smaller is —— the slower is the damping of 
2m 2m 
the oscillations. Compare these graphs with those of harmonic 


oscillations (49) shown in Fig. 12c (form (50) coincides with 


formula (49) for Pe = 0} 
2m 
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Note in addition that for large values of the factor b (for 
b> 2)/mk) formula (50) is replaced by another formula. In this 
case the weight will pass through the equilibrium position not more 
than once and after that will slowly approach this position all 
the time remaining on the same side of it (below or above). 
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Some Other Applications 
of the Concept of the Derivative 


Maximum and Minimum Values 


Let us consider some variable quantity y whose value depends 
on some other quantity x. When you say that y depends on x, 
or that y is a function of x, you mean that there is a quite 
definite value of y to correspond to every value x assumes..For 
instance, the area of a circle is a function of its radius, i.e. the 
area of the circle depends on its radius. Sine, cosine, tangent, etc. 
depend on the value of an angle, i.e. they are functions of the 
angle. These functions are termed trigonometric functions. 

So let us suppose that y is a function of x. Let us pose 
the following problem: find such a value of x that y would be 
at its maximum. Before proceeding with the solution, let us introduce 
the concept of the domain of definition of a function. We shall 
discuss this concept by using some examples. 

For the first example let us take the following function. Suppose 
V is the volume of a kilogramme of water at normal atmospheric 
pressure and at a temperature t° (Celsius). Here V depends on ¢, 
i.e. Vis a function of the quantity t. Evidently, this function is 
defined only for t’s lying inside the interval 0-100°. Indeed, at normal 
atmospheric pressure water cannot have a temperature t < 0° 
(water will turn into ice) or t > 100° (water will turn into steam). 
Hence, the function V is defined only for such values of 1 that 
satisfy the inequalities t>0 and ¢ < 100. Usually the two in- 
equalities are written together: 0 <1t < 100. Thus the function V 
is defined only for 

0<r< 100 


In other words, the domain of definition of the function V is made 
up of numbers which satisfy the condition 0<1r< 100. Such 
a domain of definition is termed a number segment because, when 
numbers are represented on the number axis, all the numbers 
which satisfy the condition 0<t< 100 fill a complete segment 
of the number axis. The numbers 0 and 100 are termed the 
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terminals or end points of the number section 0 <t < 100. All the 
other points of the section are termed its internal values or internal 
points. It is peculiar to every internal value to that its number 
segment contains both points greater and smaller than to. The terminals 
of the segment do not possess that property. 

For the second example let us consider the current i passing 
through the electric circuit (shown on the diagram in Fig. 3) 
t seconds after the circuit has been closed. In this case i is 
a function of time t. The dependence of i on tf is shown by the 
formula presented on p. 30. For what values of t is the function 
i defined? Evidently, up to the instant the circuit was closed, 
i. e. for t < 0, there had been no current passing through it, therefore 
it is reasonable to consider the current i only for t > 0. Hence, 
the domain of definition of the function i will be made up of all 
the numbers that satisfy the condition t>0. Such a domain of 
definition (it may be termed a number semi-axis) has only one end 
point t= 0, all the other points are internal. 

Finally, for the third example let us consider a function y = sin x. 
It is defined for any x, i.e. the domain of definition of this func- 
tion is the entire number axis. This region has no end points. 

There are functions whose domains of definition are very 
intricate. We, however, shall consider only such functions which 
have for their domain of definition the number segment, number 
axis or number semi-axis. 

Let us now turn again to the problem of finding the maximum 
value of a function. Can a function attain its maximum value 
at an end point of its domain of definition? Of course, it can. 
Consider, by way of example, the function V discussed above, 
which defines the volume of a kilogramme of water at normal 
atmospheric pressure and at a temperature ¢°C. Since the volume 
of water increases upon heating (from 4°C upwards), it is clear 
that the function V will have its maximum value at t = 100°, ie. 
at the end point of the domain of definition. 

In many cases differentiation enables the position of the function’s 
maximum to be determined quickly. Namely, the following 
proposition holds. 

Let y be a function of the variable x. If this function attains its 
maximum value at an internal point x = of its domain of definition, 


wave. ad 
the derivative = at this point becomes zero*. 
x 


* Subject to the condition that the derivative exists. There are such 
functions in mathematics that do not have a derivative 
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Let us prove this proposition. We shall denote the value of y 
that corresponds to the value of x (chosen from the domain 
of definition of the function) by y,. We have presumed that the 
value y, the function y assumes at x=a is its maximum 
value, i. e. 


Ya 2 Vx (51) 


for any x (chosen from the domain of definition of the function). 


The derivative 4 for x =a is determined by the relation 
x 


dy lim 22th 22 
—— = lim ————_ 
dX \ror x=a h0 h 


(52) 


Let us prove that this derivative is zero. 

We shall start making h approach zero, attributing various 
positive values to it. As long as the numerator y,,,—y, of the 
fraction under the limit sign satisfies the inequality y,,, —y, <0 
(see (51)), and h> 0, the fraction under the limit sign will itself 
be non-positive (i. e. it may be zero, or negative). But in that case 
the limit of this fraction, too, cannot be positive, i. e. derivative 
(52) cannot be a positive number. 

Next let us make h approach zero, attributing various negative 
values to it. In this case, as before, y,., — y, <0 (see (51)), but 
h <0, and for this reason the fraction under the limit sign will 
itself be non-negative. But in that case the limit of this fraction 
(i. e. the derivative we are interested in) cannot be negative. 


cs hae d . 
Thus, the value of the derivative = for x =a can be neither 


positive nor negative, therefore, it must be zero, and this proves 
our proposition. *) 

In this proof the essential fact was that a is an internal point 
of the domain of definition of the function. Indeed, we attributed 
positive and negative values to h so that a+h assumed values 
greater and less than a. 

Suppose a is an end point. Then the domain of definition 
includes either values greater than a, or less than a, i. e. the above 
proof is inapplicable. 

When dealing with the problem of the minimum (and not maximum) 
value of a function, the analysis is quite similar. As a result we 


* It would certainly not be out of place to point out that if the limit 
exists, its value must be the same for h>O and h<O, ie. 
lim = lim — translator's note. 
h>O h<O 
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shall prove that if a function assumes its minimum value at an 
internal point of its domain of definition, the function’s derivative 
at this point will become zero. Combining both these cases — 
of maximum and minimum values — we obtain the Fermat’s theorem, 
named after a French mathematician of the XVIIth century. 

THEOREM If a function assumes its maximum (or minimum) value 
at an internal point of its domain of definition, the function's 
derivative at this point becomes zero. 

This theorem serves as a basis for obtaining maximum and 
minimum values with the aid of differentiation. We must find the 
derivative of the function in question and the corresponding 
internal points of its domain of definition where the derivative 
becomes zero. One should look for the point where the function 
attains its maximum (or minimum) value among these points 
(where the derivative becomes zero) or among the end points of 
the domain of definition. 

Example 11. A battery having an electromotive force E and 
an internal resistance r is connected to the terminals of a conductor 
(for instance, a heater). What should the conductor’s resistance be 
so that it draws the maximum power from the battery? 

Solution. Denote the conductor’s resistance by R. In this case 
the total resistance of the circuit will be R + r, and, accordingly, 


the current passing through the circuit will be i = — 
The power supplied by the battery to the conductor is expressed 
2 
by the formula W = i?R, i.e. W= _ ER _ 53 
y i7R, ie Rtn? (53) 


Hence, problem may be formulated as follows: for what R does 
function Was expressed by formula (53) attain its maximum value? 

The domain of definition of the function W is the semi-axis 
R20 (we consider only the case of positive ohmic resistance). 


Find the derivative 4. 
aR 


E(R+h) ___E?R 
dW _iim (R+eh+rP (R+rP 
“dR hd h ee 
— tim EZUR +) (R+ rh - R(R+ AHP) _ 
h—0 h(R+h+ryP(R+r/) 
‘ r? — R?—Rh r? — R? r—R 
= lim 


no(R+h+rP(R+rP (R+r (Rr 
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lor the derivative ss i.e. the fraction 
r-—R 
(R +r) 


to become zero, its numerator r—R should become zero, i.e. 
it should be R =r. 

Hence the power W can attain its maximum (or minimum) value 
either for R=r or at the end point of its domain of definition 
R=0. For R =0 the power W is zero, too (this is the minimum 
value). Therefore the power can attain its maximum value only 
for R=r, i.e. when the conductor's resistance is equal to the 
internal resistance of the battery. 

Shall the power really be at its maximum for R =r? Indeed, 
we have proved only that the power may attain its maximum 
value for R=r, but this in itself does not mean that it really 
will be so. 

It is not difficult to make sure that when R =r, the power W 
is really at its maximum. Indeed, when R =0, the power W is 
zero, too; when R is very great, the current i will be very smail 
and hence the power will be small (for the voltage drop across 
the conductor’s terminals does not exceed £). It is therefore clear 
that the power must attain its maximum value for some (not 
very great) value of R. But since the power must attain its maximum 
value (and it may do so only for R=r), therefore it stands 
to reason that we do, indeed, obtain the maximum power when 
R=r. 

Example 12. The problem is how to build a steam _ boiler 
in the shape of a cylinder so that it will have the required 
volume V. It is desirable to keep the total surface of the boiler 
down to a minimum (in that case the minimum amount of metal 
will be- required to build the boiler; besides the smaller is the 
surface of the boiler the less are the heat losses due to contact 
with the atmosphere). Find the best (optimum) dimensions of the 
boiler. 

Solution. Denote the radius of the cylinder’s base by R and its 
height by h. Then 


V =nRh 
i.e. , 
h = 2 


The cylindrical surface is S = 2nR? + 2aRh, i.e. 
Saar ae (54) 
R 
We must find out for what values of R the quantity S (which 
depends on R, i.e. is a function of radius R) assumes its minimum 


value. Find the derivative aS 
dR 


2V 2V 
2 = 2, 4 
[2ne +m haere Ee + x 


aS 


= li 
roar h 
2Vh 
5) aed AOD 
2n (2Rh + h*) R(R+h 
= lim = 
4-0 h 
: 2Vv 2V 
= im | 4nR + 2nh — Rein | = 4nR — RE 
3 
Equating the derivative AS to zero, we find: R = Be and 
dR : 2n 
therefore 


3. —— 
V 4V V 
h= pr = iio er 

In other words, the height of the cylinder should be equal to its 
diameter. 

Do we really obtain the minimum value of the cylindrical 
surface for that value of R? This we can easily check. Indeed, 
for very great values of R the surface S is also very great (since 
the value of the first term in the expression for S is great — see 
(54)). For very small values of R the value of the surface S is 
also very great (this time it is the second term that accounts for 
this). Therefore for some intermediate value (not too great and 
not too small) of R the quantity S must assume its minimum value. 


But since the derivative = becomes zero only for one value of 


R, this is the value of R to which the minimum surface of the 
cylinder corresponds. 

We shall confine ourselves to these two examples. If instead, 
the reader will find many such examples in textbooks and problem 
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books. The reader may be advised to solve some of them provided 
he does not ignore the final stage, i.e. the proof that there is 
indeed a maximum or a minimum value at the point obtained. 
More perfect methods are described in higher mathematics courses 
which enable the decision to be made whether the function does 
indeed assume its minimum or maximum value at the point 
obtained. Besides, there are rules for calculating derivatives. 
The author worked on the assumption that the reader is not 
familiar with these rules, and for this reason the derivatives 
in the examples above have been obtained with the aid of direct 
calculations. 


The Problem of Drawing a Tangent 


Let L be a curve and Mo a point on it. Let us discuss the 
problem of drawing a tangent to the curve L at point Mo. First 
let us discuss how the tangent is defined in mathematics. Choose 
point M also lying on the curve Land draw a straight line MyM 
which we shall term a secant, for it intersects the curve L at least 
at two points My and M. Should point M move along the curve 
Lapproaching point My (Fig. 13 shows the sequence of positions 


Fig 13 


M, M’', M", ... of point M), the secant M)M would turn around 
point Mo. If the secant MoM in turning tends to coincide with 
some straight line M)K, this limiting line M)K is termed the 
tangent to the curve L at point Mo. 

Suppose now that the curve L is drawn in a plane containing 
a coordinate system so that each point M of the curve L has its 
corresponding x- and y-coordinates. Denote the x-coordinate of 
point M, by a (Fig. 14), and the length of section NoN by h. 
Then the x-coordinate of point M will be a+h. We denote 
the y-coordinate of point My by y,, and that of point M by 
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Fig 14 
Ya+h- The length of section MP is 
MP = MN — PN = MN — MWNo = Yash — Ya 


and we therefore have 
MP MP as Vat+h — Ya 
MoP  NoN | h 


Denote an angle PMoK, i.e. the angle between the x-axis and 
the tangent, by «. Then as M approaches Mo, i.e. as a section 
NoN =h tends to zero, the angle PM )M will approach a, and 
the tangent of the angle PM,)M tana. Hence, we obtain from 
relation (55) in the limit (for h— 0) 


tana = lim Yarn Yo _ dy 
4-0 h dx for x=a 


tan <2 PM)M = (55) 


It follows that the slope of the tangent line is equal to the value 
of the derivative of the y-coordinate with respect to the x-coordinate 
for x =a, where a is the point of tangency. 

Example 13. Consider a sine curve (Fig. 15), i.e. a curve 
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whose x- and y-coordinates are related by the equation 
y=sinx 


How should a tangent line be drawn to this curve at some point 
Moy whose x-coordinate is equal to a? 
We know already how to find the slope of this tangent line 
dy 


d . 
tana = —— =—— sinx = cos a 
dx for x=a dx 


for x=a 


(see (36)). Hence, to draw a tangent line we must find cosa 
(which is quite easy since we know the length of section 
MoNo =y=Sina) and draw a straight line M)K so that 
tan « = cos a. For instance, when a = 0, we obtain: tana = cos 0 = 1, 
i.e. the tangent line to a sine curve drawn through the origin of 


- rt ‘ : rv 
coordinates makes an angle a with the x-axis. For t= We 


: tT 1 : ; ; 
obtain: tan « = cos san whence, using the trigonometric tables, 


we obtain that the angle «, which the tangent line makes with 
the x-axis, is equal approximately to 26°34’. 


Modelling 


We have seen above that various physical phenomena which 
at first glance have nothing in common may be described by 
the same differential equations. That was the case with a body 
falling in a resistant medium, with the current being switched 
on in a circuit and with radioactive decay. This was the case with 
the three problems leading to harmonic oscillations. But when two 
phenomena are described by one equation, then the solution of 
this equation describes both phenomena. In other words, both 
phenomena take the same course (i.e. the values of the corres- 
ponding physical quantities change identically with time). Therefore 
we may, for example, study the electromagnetic oscillations 
in a circuit by observing the oscillations of a pendulum, and vice 
versa. This very simple remark proves to be very significant. 

Suppose we have a design of an intricate machine before us 
and want to check whether the design was calculated correctly. 
It takes much time and money to build the machine, and still 
we must make sure beforehand that the design is the right one. 
Let us write out the equations of motion for the machine (these 
are usually differential equations). Should we solve these equations 
we would know how the machine is going to work and whether 
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it has been correctly designed. But it often proves easier to act 
as follows. Let us build another device (for instance, assemble 
an electric circuit) which would be described by the same differen- 
tial equations as the machine to be designed. If we are able to do 
it, it will be sufficient to study the behaviour of the device we 
have built to learn how the machine we are designing is going 
to work. Thus, this device is a model of machine. 

Such “modelling”, too, uses the concept of the derivative as its 
basis, for to be able to build the right model one must know 


the differential equations which describe the operation of the 
machine. 


Afterword 


The concepts of the derivative and the differential equation are 
used extremely wide in mathematics, physics, astronomy, and 
engineering. The discipline that studies the properties and the 
applications of these concepts is “higher mathematics”. Unfortu- 
nately, within the space of this booklet it would have been 
difficult to explain ideas underlying the definition of integration® 
which, in a sense, is an operation reciprocal to differentiation and 
together with the latter forms the cornerstone of higher mathe- 
matics. 

The author would like to draw the reader’s attention to the 
point that the concepts of higher mathematics, including the 
concept of the derivative, are not abstract, but constitute 
a mathematical reflection of processes taking place in nature (for 
instance, of the speed of mechanical motion). These concepts 
developed in close contact with problems posed by everyday 
life — primarily in the fields of mechanics (the problem of calculating 
velocities) and geometry (the problem of drawing a_ tangent). 
F. Engels wrote: “Mathematics together with other sciences arose 
out of man’s requirements.” This statement fully applies to higher 
mathematics. The concept of the derivative created in conjunction 
with the study of the motion of bodies and the variation of 
quantities itself reflects this motion, for it treats a variable quantity. 
“The turning point in mathematics was Descartes’ variable quantity. 
Thanks to it mathematics was invaded by motion and dialectics 
and it, too, created an immediate need for differential and integral 
calculus the rudiments of which were soon to be established and 
which was completed as an entity but not discovered by Newton 
and Leibnitz” (Engels). 


*) See I. P. Natanson Summirovanie beskonechno malykh velichin (The 
Summation of Infinitesimals) published in this series (No. 12). 
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Thus the concepts of higher mathematics arose out of man’s 
requirements, primarily in connection with the study of mechanical 
motion. But motion in nature is not only limited to its mechanical 
form. Time and again new phenomena reveal themselves to 
scientists, and mathematics has to describe newly discovered 
phenomena and forms of motion. Physical theories created in recent 
times (the theory of relativity, quantum physics, the theory of the 
nucleus) require new mathematical apparatus. Some mathematical 
disciplines were born during the last decades. Mathematics is not 
an abstract science but a science closely connected with life which 
develops together with the growth of our knowledge of the physical 
world. 
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